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^ ' Abstract 

o ■ 

O ' Multi-electron production is studied at high electron transverse momentum in positron- 

and electron-proton collisions using the HI detector at HERA. The data correspond to an 
_ integrated luminosity of 115 pb^^. Di-electron and tri-electron event yields are measured. 

■ Cross sections are derived in a restricted phase space region dominated by photon-photon 
collisions. In general good agreement is found with the Standard Model predictions. How- 

g_) I ever, for electron pair invariant masses above 100 GeV, three di-electron events and three 

■ tri-electron events are observed, compared to Standard Model expectations of 0.30 ± 0.04 
and 0.23 ± 0.04, respectively. 
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1 Introduction 



In this paper we describe the first measurement of multi-electron production at high transverse 
momentum (P^) in electron ^-proton interactions at HERA. Within the Standard Model (SM), 
the production of multi-lepton events in ep collisions proceeds mainly through photon-photon 
interactions; photons radiated from the incident electron and proton interact to produce a pair of 
leptons, 77 Q. At large invariant masses, multi-lepton production may be sensitive to 

new phenomena, for instance the production of a doubly charged Higgs boson or processes 
involving bileptons, generic bosons carrying two units of lepton number lISj. 

The analysis presented here is based on data recorded in 1994-2000 by the HI experiment. 
Electrons of 27.6 GeV collided with protons of 820 or 920 GeV, corresponding to centre-of- 
mass energies y/s of 301 GeV or 319 GeV, respectively. The total integrated luminosity of 
115.2 pb^^ consists of 36.5 pb~^ of e~^p collisions recorded at ^/s = 301 GeV and 65.1 pb^^ 
recorded at 319 GeV, as well as 13.6 pb^^ of e^p collisions recorded at = 319 GeV. A related 
study of muon pair production is presented in 



2 Standard Model Processes and their Simulation 

The main SM processes contributing to multi-electron production at HERA are summarised in 
figure [TJ The dominant contribution, shown in diagram [1^, is due to electron pair production 
through the interaction of two photons radiated from the incident electron and proton. Electron 
pairs can also originate from internal conversion of a photon (7) or a Z° boson, radiated either 
from the electron line (diagram [1J5) or from the quark line (diagram [it). The pole due to the 
electron propagator in diagrams [it and[lj5 corresponds to an e+e^ interaction in which one 
of the electrons is emitted from a photon radiated from the proton. This mechanism is called 
the Cabibbo-Parisi process. Its contribution is one order of magnitude lower than the photon- 
photon contribution, except at high transverse momentum, where it is more significant due to s- 
channel boson production (diagram[Tj5). In diagram[lt, the pole due to the quark propagator 
corresponds to the Drell-Yan process, qq e+e^. Its contribution is small compared with the 
photon-photon and Cabibbo-Parisi processes ||5J. 

If the photon coupled to the incoming electron has a high virtuality, the incident electron 
can scatter through a large angle and with high transverse momentum. If the scattered electron 
is observed in the detector, it is indistinguishable from the pair-produced electron of the same 
charge. 

The Monte Carlo generator GRAPE |6| simulates lepton pair production in ep collisions 
using the full set of electroweak matrix elements at the tree level, with the exception of the 
Drell-Yan pole contribution. GRAPE is based on the automatic Feynman graphs calculation 
program GRACE |7 1. Initial and final state radiation processes (QED and QCD parton showers) 
are simulated in the leading log approximation. The production of r lepton pairs and their 
subsequent electronic decay is also simulated with GRAPE and composes about 2% of the 
multi-electron event sample. The Drell-Yan contribution is simulated using the PYTHIA JSJ 

'in this paper the term "electron" is used generically to refer to both electrons and positrons. 
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event generator. Its contribution is found to be very small in the phase space relevant here and 
is neglected in the following. 

In GRAPE the proton interaction is divided into in three phase space regions: elastic, quasi- 
elastic and inelastic. In the case of elastic scattering, ep ee^e^p, the proton vertex is de- 
scribed in terms of dipole form factors. The quasi-elastic domain is defined by requiring that 
the mass of the hadronic final state Mx < 5 GeV or that the virtuality of the photon coupled to 
the proton Qp<l GeV^. In the region Mx < 2 GeV a resonance parameterisation is used 
for the proton vertex. In the remaining quasi-elastic phase space, a fit to photoproduction and 
deep inelastic scattering data is used [10|. The inelastic regime corresponds to electron-quark 
interactions with Mx > 5 GeV and Q"^ > 1 GeV^. In this case, the proton structure is pa- 
rameterised using the CTEQ5L parton distributions ifTTIl . The fragmentation and hadronisation 
processes are simulated using the SOPHIA program [12J in the quasi-elastic and PYTHIA lIHl 
in the inelastic regime. 

The uncertainty attributed to the GRAPE calculation in this analysis is 3%. This value 
results mainly from the uncertainties in the QED matrix element calculation (1%), the 
parameterisation of the structure functions and the phase space separation between quasi-elastic 
and inelastic processes. 

The GRAPE prediction is cross-checked using the LPAIR generator |[T3, which contains 
only the photon-photon process. When restricted to this process, the total and differential 
cross sections produced with LPAIR and GRAPE agree at the percent level. The additional 
diagrams in GRAPE increase the predicted cross section by 10 % on average in the phase space 
considered here. The increase is more pronounced (up to 30 %) for di-lepton masses which are 
either very low (photon internal conversions) or around 90 GeV (Z" resonance production). 

The main experimental backgrounds to multi-electron production are processes in which, 
in addition to a true electron, one or more fake electrons are reconstructed in the final state. 
The dominant contribution arises from neutral current Deep Inelastic Scattering (DIS) events 
{ep — > eX) in which, in addition to the scattered electron, hadrons or radiated photons are 
incorrectly identified as electrons. QED Compton scattering {ep e'jX) can also contribute if 
the photon is misidentified as an electron. The DIS and elastic Compton processes are simulated 
using the DJANGO lHUl and WABGEN ^T5^ generators, respectively. 

All generated events are passed through the full GEANT ifT^ based simulation of the HI 
apparatus and are reconstructed using the same program chain as for the data. 

3 Experimental Conditions 

A detailed description of the HI detector can be found in ifTTll . The components essential for 
this analysis are described briefly here. 

A tracking system consisting of central and forward^ drift chambers is used to measure 
charged particle trajectories and to determine the interaction vertex. The central tracker is com- 
posed of two concentric cylindrical drift chambers with an active detection region starting at a 

^The origin of the HI coordinate system is the nominal ep interaction point. The direction of the proton beam 
defines the positive z-axis (forward direction). Transverse momenta are measured in the x-y plane. Polar (6) 
and azimuthal (0) angles are measured with respect to this reference system. The pseudorapidity is defined as 
77 = -logtaii(6'/2). 
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radius of 22 cm. The angular range 37° < 6 < 144° is covered by both chambers. The inner 
drift chamber provides full acceptance for particles in the range 22° < 9 < 159°. Particles at 
9 = 20° cross 83% of its acceptance region. Transverse momenta (Pt) are determined from 
the curvature of the particle trajectories in a magnetic field of 1.15 Tesla. The central tracking 
system provides transverse momentum measurements with a resolution of ap^/ = 5 x 10^^ 
GeV~^. The tracking is complemented in the region 7° < 9 < 25° by a system of drift chambers 
perpendicular to the beam axis. 

Hadronic and electromagnetic final state particles are absorbed in a highly segmented liquid 
argon calorimeter lITHI covering the range 4° < 6* < 153°. The calorimeter is 5 to 8 hadronic in- 
teraction lengths deep, depending on the polar angle and has an electromagnetic section which 
is 20 to 30 radiation lengths deep. Electromagnetic shower energies are measured with a pre- 
cision of a{E)/E = 12%/^E/GeV © 1% and hadronic shower energies with a precision of 
a{E)/E = 50%/ y^E/GeV © 2%, as measured in test beams lfT9ll . The electromagnetic energy 
scale is known to 0.7% in the central region and to 3% in the forward region. The hadronic 
energy scale is known to 2%. 

In the backward region, energy measurements are provided by a lead/scintillating-fibre 
calorimeter^ |i20| covering the range 155° < 9 < 178°. The calorimeter system is surrounded 
by a superconducting coil with an iron yoke which is instrumented with streamer tubes. The 
electron and photon taggers located downstream of the interaction point in the electron beam 
direction are used to determine the luminosity through the measurement of the Bethe-Heitler 
ep e'jp process. 

The trigger used relies on the liquid argon calorimeter signals and has an efficiency which 
is greater than 95% for events in which an electron of energy above 10 GeV is detected. 



4 Data Analysis 

4.1 Multi-electron event selection 

The multi-electron event selection is based on a procedure which is designed to minimise the 
contribution of fake electrons, while keeping a high efficiency for identifying true electrons and 
allowing reliable monitoring of the overall selection efficiency. 

As a first step, electron candidates with energies above 5 GeV are identified in the liquid 
argon and backward calorimeters, in the range 5° < < 175°. Electromagnetic showers are 
identified with an efficiency of better than 98% using pattern recognition algorithms based on 
the geometric profiles expected for electrons. The remaining calorimeter clusters are attributed 
to hadronic activity and are combined into jets using an inclusive k-r algorithm lETl . with a 
minimum jet transverse momentum of 4 GeV. Electron candidates are required to be isolated by 
demanding that they are separated from other electrons or jets by at least 0.5 units in the t] — (p 
plane. In addition, the total hadronic energy within 0.75 units in — of the electron direction 
is required to be below 2.5% of the electron energy. 

^This device was installed in 1995, replacing a lead-scintillator "sandwich" calorimeter 1171 . 
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In the region of angular overlap between the liquid argon calorimeter and the central drift 

chambers (20° < 9 < 150°), the calorimetric electron identification is complemented by track- 
ing conditions. In this region it is required that a high quality track be geometrically matched to 
the electromagnetic cluster with a distance of closest approach to the cluster centre of gravity 
of less than 12 cm. The starting radius of the measured track, defined as the distance between 
the first measured point in the central drift chambers and the beam axis, is required to be be- 
low 30 cm in order to reject photons that convert in the central tracker material beyond this 
radius. The transverse momentum of the associated track P^,*'' and the calorimetric transverse 
momentum P|. are required to satisfy the condition 1/P|," — 1/P|. < 0.02 GeV~^. No other 
high quality track is allowed within 0.5 units in — of the electron direction. These addi- 
tional constraints strongly reduce the contribution of fake electrons from misidentified photons 
and hadrons. The resulting electron finding efficiency is 88%. Electrons selected in this polar 
angular range are called "central electrons" hereinafter. 

Due to the higher material density in the forward region (5° < ^ < 20°) the electrons are 
more likely to shower and therefore no track conditions are required. The same applies in the 
backward region (150° < 9 < 175°). The forward electron energy threshold is raised to 10 GeV 
in order to reduce the number of fake electrons arising from hadrons in DIS events. 

The final multi-electron event selection requires that there be two central electron candi- 
dates, of which one must have P|. > 10 GeV and the other P^ > 5 GeV. Additional electron 
candidates are identified in the central and backward regions with > 5 GeV and in the for- 
ward region with > 10 GeV. The electron candidates are ordered according to decreasing 

The selected events are classified as "2e" if only the two central electron candidates are 
identified and "3e" if exactly one additional electron candidate is identified. A subsample of 
the "2e" sample, labelled "77", is selected in order to measure the pair production cross section 
in a well defined phase space region dominated by photon-photon collisions with low back- 
ground. In this subsample, the two electrons must be of opposite charge and a significant deficit 
compared to the initial state must be observed in the difference E — P^of the energy and longi- 
tudinal momentum of all visible particles (E — P^ < 45 GeV) These two conditions ensure 
that the incident electron is lost in the beam pipe after radiating a quasi-real photon of squared 
four-momentum Q"^ lower than 1 GeV^. 

4.2 Background studies 

DIS and Compton processes can contribute to the selected multi-electron sample if a photon or 
a particle from the hadronic final state is misidentified as an electron. In order to quantify the 
uncertainty on the background prediction and test the performance and reliability of the electron 
identification procedure, several samples in which these background processes are enhanced are 

studied. 

^For fully contained events or events where only longitudinal momentum along the proton direction (+z) is 
undetected, one expects E — = 2E^ = 55.2 GeV, where E° is the energy of the incident electron. If the 
scattered electron is undetected, the threshold E — < 45 GeV corresponds to a cut on the fractional energy loss 
y={E- P,)/2El < 0.82. 
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• Electron misidentification in the central region is investigated by measuring the probabil- 
ity of selecting a second electron in addition to the scattered electron in DIS candidate 
events when the track quality criteria are relaxed. From a sample of 244980 DIS candi- 
dates, 1563 events with a second electromagnetic central cluster are selected if no tracking 
conditions are applied to this cluster. The fake electrons are predominantly photons from 
DIS or Compton events (figure|2ti). This contribution is greatly reduced by requiring a ge- 
ometrical track-cluster match (figure |2b)- The remaining sample consists of 250 events. 
The fake electron background is described by the simulation at the 20% level. 

• Electron misidentification in the forward region affects only the "3e" selection. The dom- 
inant contribution to this background is the misidentification of a hadron as an electron. 
The fake electron background in this region is studied with the DIS event sample in an 
analysis similar to that described above. Fake electrons in the forward region are also 
searched for in an inelastic Compton event sample with one electron and one photon in 
the central region. The fake electron background in the forward region is described by 
the simulation at the 20% level. 

• Detailed studies of photon conversions are performed using a sample enriched with elastic 
Compton events, selected by requiring one central electron plus a second central electro- 
magnetic cluster (photon candidate) and no significant additional energy in the calorime- 
ters. Distributions of the charged tracks associated with the photon candidate are shown 
in figures |2t and |2|l. The number of tracks (figure |2l:) and their starting radius (fig- 
ure |2|l) are well reproduced by the simulation. In figure |2|l, the central tracker structure 
is visible as peaks in the distribution, corresponding to photon conversions in the tracker 
walls. The first peak is populated by tracks associated with true electrons and by tracks 
from conversions which occur before the active tracker volume. The second peak is due 
to photon conversions in the dead material between the inner and outer central trackers. 
These conversions are described by the simulation to better than 20%. 

Based on those studies, the uncertainty on the background simulation is estimated to be 20%. 

4.3 Systematic uncertainties 

The systematic uncertainties are related to the measurement of the electron pair production 
process, to the theoretical description of this process and to the background simulation. 

The main measurement uncertainty is due to the tracking conditions in the electron iden- 
tification procedure. The electron track association efficiency is measured with a DIS sample 
selected with E — > Ah GeV and a single electromagnetic cluster in the calorimeter with a 
transverse momentum above 10 GeV in the polar angle interval 20° < 6 < 150°. The measured 
average track association efficiency is 90% and varies only slightly with the track momen- 
tum and polar angle. This efficiency, measured to a precision ranging from 3% for polar angles 
around 90° to 15% at the forward edge of the angular acceptance of the central tracker {6 = 20°), 
is well described by the simulation. Uncertainties on the energy scales of the calorimeters, on 
the trigger efficiency and on the luminosity measurement are also taken into account. The total 
measurement uncertainty is typically 7% for the "2e" selection and 10% for the "3e" selection. 
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The theoretical uncertainty on the pair production process cross section, calculated with 
GRAPE, is 3%, as explained in sectional The uncertainty on the Compton and DIS background 
contributions is 20%, deduced from the studies presented in section 

The error on the event yields predicted by GRAPE (section 1511) contains all measurement 
and theoretical errors added in quadrature. The uncertainty on the total SM prediction also 
includes the errors on the Compton and DIS backgrounds. 

The error on the extracted cross sections (section l5?2b includes all measurement and back- 
ground errors as described above. Theoretical errors are applied to the GRAPE prediction of 
the cross sections. 



5 Results 

5.1 Multi-electron event samples 

The multi-electron event yields are summarised in table [H The observed numbers of "2e" and 
"3e" events are in agreement with the expectations, as is the number of events in the "77" 
sample. No event is found with four or more identified electrons. 

The distributions of longitudinal momentum balance E — P^, missing transverse momentum 
pmtss hadronic transverse momentum p^i^-d-rons presented in figure |3l The "3e" events 
accumulate at E — values around 55 GeV, as expected if the scattered electron is visible in 
the detector. The "2e" events show a tail at lower E — P^, due to the scattered electron being 
lost in the beam pipe, corresponding to the dominant 77 topology. The missing transverse 
momentum p™*'^* is taken to be the modulus of the vector sum of the transverse momenta of 
all visible particles. The P™'^'* distributions are consistent with the expectation for no emission 
of undetected particles with substantial transverse momentum. The spectrum of the transverse 
momentum p}j;°-'^'^°^^ of all visible particles except identified electrons is also well described by 
the SM prediction. 

The distributions of the individual electron transverse momenta P|.' are steeply falling as 
shown in figure IH The "2e" and "3e" samples are in good overall agreement with the SM, 
except for three "2e" events with P^^ above 50 GeV, where the SM expectation is small. 

The distribution of the invariant mass of the two highest Pr electrons in the event (M12) 
and the correlation with the scalar sum P^^ + P^ are shown in figure 13 The agreement with 
the SM prediction is good at low M12. However, three "2e" and three "3e" events are seen 
with invariant masses M12 above 100 GeV, where the SM expectation is small. The three "2e" 
events are the same as those observed at high P^\ The invariant masses Afi3, M23 and M123 
of the other possible electron combinations in the "3e" sample are shown in figure |6l No event 
is seen with an unexpectedly high M13 or M23. The three high M12 events also give rise to 
the largest tri-electron masses M123. The comparison of the observed events with masses M12 
above 100 GeV with the SM expectations is presented in table|2| These events are discussed in 
detail in section 1531 
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5.2 Cross section measurement 

Using the selected "77" sample, electron pair production cross sections are measured in the 
kinematic region defined by 

20° < r^'^ < 150°, P^' > 10 GeV, P^^ > 5 GeV, y < 0.82, < 1 GeV^ . 

For this measurement, the data samples collected at ^/s=30l GeV and 319 GeV are combined 
taking into account their respective luminosities. Assuming a linear dependence of the cross 
section on the proton beam energy, as predicted by the SM, the resulting cross section corre- 
sponds to an effective ^/s = 313 GeV. 

After background subtraction, the observed number of events is corrected for acceptance 
and detector effects to obtain the cross section. The generator GRAPE is used to calculate the 
detector acceptance A for this region of phase space. The acceptance accounts for detection 
efficiencies and migrations. The cross section is thus 

Ndata — Nbgd 

CA ' 

where N^ata is the number of events observed, Nbgd is the number of events expected from 
background processes (DIS and Compton) and C is the integrated luminosity of the data sample. 

The ep — > ee'^e'X cross section, integrated over the phase space defined above, is 

a = (0.59 ±0.08 ±0.05) pb, 

where the first error is statistical and the second systematic, obtained as described in section 4.2. 
This result agrees well with the SM expectation of (0.62 ± 0.02) pb, calculated with GRAPE. 
The differential cross sections as a function of M12 and p^^'^^o"* are shown in figure|7]and 
table |3l The measurements are in good agreement with the expected cross sections. 

5.3 Discussion of high mass events 

All six events with M12 > 100 GeV were recorded during positron-proton collisions. For these 
events, displayed in figureEl all available detector information supports the interpretation of the 
electron candidates as being true electrons. The electromagnetic shower shapes are checked in- 
dividually and found to be similar to those expected from the calorimeter response to electrons. 
All central tracks yield a specific ionisation in the central drift chamber as expected for single 
electrons. The measurements of the central electron momenta by the tracker and the calorimeter 
are compatible within errors. The forward electron candidates in the "3e" events 4, 5 and 6 all 
have at least one track pointing to the calorimetric energy cluster, although no such requirement 
is made in the identification procedure. 

Although classified as "2e", event 1 also contains a third electron candidate with energy 
below the identification threshold. Similarly, event 3 has a compact electromagnetic energy 
deposit located at the forward boundary of the liquid argon calorimeter, outside the electron 
identification fiducial volume. This event also contains a low energy converted photon radiated 
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close to electron e^. With the exception of event 6, which has a high energy forward hadronic 
jet, the events show no hadronic activity in the detector. It should be noted that the high mass 
di-electron topology differs for the observed events classified as "2e" and "3e". In the former 
case, the high-mass is formed from two central high-P^ electrons, whereas in the latter it is 
formed from one forward and one central electron, both of intermediate Pt (figure |5l). 

The event kinematics of the six high mass events are presented in table |4| The electron 
energy, electron polar angle and forward and backward electron azimuthal angles are mea- 
sured from the calorimetric deposits. For the central electrons, the azimuthal angle is measured 
from the associated track, which yields a better precision. The electric charge of the electrons, 
measured in the central region from the associated track curvature, is given in tablelHif the sig- 
nificance of its determination exceeds two standard deviations. All events are compatible with 
the presence of one and two in the final state, as expected from pair production process. 

Imposing longitudinal and transverse momentum conservation, a constrained fit can be per- 
formed to improve the kinematic measurement. This corresponds to an adjustment of particle 
observables (energy, polar and azimuthal angles) within experimental errors in order to achieve 
E — Pz = 55.2 GeV and p™*'^^ = GeV. In event 2 only two electrons are visible and the 
measured E — P^ value is significantly lower than 55 GeV. For this event it is supposed that 
the scattered electron has escaped down the beampipe and therefore the E — P^ constraint is 
removed in the kinematic fit. The Mu values obtained from the kinematic fit are indicated in ta- 
blelH They are consistent with the direct measurements. The errors are reduced by more than a 
factor of two with the exception of event 2. The per degree of freedom of the kinematic fit is 
in the range 0.7 to 1.7 for the six events, showing that the kinematics of the six high mass events 
are well understood within the measurement errors. The Mu values are incompatible with the 
interpretation of the six high mass electron pairs as the decay of a single narrow resonance. The 
same is true for the Mus values in the "3e" high mass events. 

6 Summary 

High-P^- multi-electron production is measured for the first time in ep scattering at HERA. 
The di-electron and tri-electron event yields are in good overall agreement with the SM pre- 
dictions. No events are observed with more than three identified electrons, again in agreement 
with the SM expectation. Differential cross sections for electron pair production are extracted 
in a restricted phase space region dominated by photon-photon interactions and are found to 
agree with the predictions. 

Within the di- and tri-electron samples, the invariant mass of the two electrons with the 
highest transverse momenta is studied. For masses above 100 GeV, three events classified as di- 
electrons and three events classified as tri-electrons are observed, compared to SM expectations 
of 0.30 ± 0.04 and 0.23 ± 0.04, respectively. 
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Selection 


Data 


SM 


Pair Production (GRAPE) 


DIS + Compton 


"2e" 


108 


117.1 ±8.6 


91.4 ±6.9 


25.7 ±5.2 


"3e" 


17 


20.3 ±2.1 


20.2 ±2.1 


0.1 ±0.1 


"4e" or more 





0.12 ±0.04 


0.12 ±0.04 


< 0.02 (95% C.L.) 


"77" subsample 


42 


44.9 ±4.2 


43.7 ±4.2 


1.2 ±0.4 



Table 1 : Observed and predicted multi-electron event yields for the samples described in the 
text. The analysed data sample corresponds to an integrated luminosity of 115.2 pb~^ The 
errors on the predictions include model uncertainties and experimental systematic errors added 
in quadrature. 



Selection 


Data 


SM 


Pair Production (GRAPE) 


DIS + Compton 


"2e" M12 > 100 GeV 
"3e" M12 > 100 GeV 


3 
3 


0.30 ±0.04 
0.23 ±0.04 


0.21 ±0.03 
0.23 ±0.03 


0.09 ±0.02 
< 0.02 (95% C.L.) 



Table 2: Observed and predicted multi-electron event yields for masses M12 > 100 GeV for the 
samples described in the text. The analysed data sample corresponds to an integrated luminosity 
of 115.2 pb~^. The errors on the predictions include model uncertainties and experimental 
systematic errors added in quadrature. 
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Variable range 
[GeV] 


Measured cross section 
[pb/GeV] 


Pair Production (GRAPE) cross section 
[pb/GeV] 




10 < < 15 


0.092 ± 0.016 ± 0.009 


0.090 ± 0.003 


15 < P^' < 20 


0.021 ± 0.008 ± 0.002 


0.023 ±0.001 


20 < P^' < 25 


0.0053 ±0.0037 ±0.0007 


0.0065 ± 0.0002 


da/dMi2 


15 < Mi2 < 25 


0.030 ± 0.007 ± 0.003 


0.027 ±0.001 


25 < Mi2 < 40 


0.015 ± 0.004 ± 0.001 


0.017 ±0.001 


40 < Mi2 < 60 


0.0020 ±0.0012 ±0.0002 


0.0026 ±0.0001 


da / dP!^"^'''^' 


< P!^<'drons ^ 5 


0.079 ± 0.014 ± 0.009 


0.087 ±0.003 


5 < P^d.^<^^ < 12 


0.028 ±0.011 ±0.002 


0.018 ±0.001 


12 < p^«<^'-o"« < 25 


0.0032 ±0.0023 ±0.0005 


0.0041 ± 0.0001 



Table 3 : Differential cross sections da / dP!j} , da / dMi2 and da / (ip^'^*""'' for the process ep 
ee^e^X in a restricted phase space (see text). The differential cross sections are averaged over 
the quoted intervals. The first error is statistical and the second is systematic. Theoretical 
predictions with GRAPE are also shown. 





P X P X 

a) b) 




Figure 1 : The main processes involved in lepton pair production. Examples of Feynman di- 
agrams are shown for: a) photon-photon interaction; b) and c) 'j/Z'^ boson conversion. The 
hadronic final state (X) can be a proton (elastic process) or a higher mass system (quasi-elastic 
and inelastic processes). 
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Particle 


E [GeV] 


6 [degrees] 


(jf) [degrees] 


Charge (significance) 


Run 83507 
Event 16817 


Multi-electron Event 1 (2e) 

£;-P^=54.0± 1.1 GeV p™'"'=3.1 ± 1.8 GeV 
Mi2 = 1 1 1 .2 ± 2.4 GeV = 1 1 1 .3 ± 0.4 GeV 


ei 
e2 

low energy e 


90.3 ±3.1 
53.6 ± 1.4 
4.4 ± 0.3 


36.6 ± 0.2 
69.6 ± 0.3 
44.3 ± 0.3 


98.48 ± 0.05 
-77.05 ± 0.05 
-155.46 ±0.03 


-(4a) 
undetermined 
+ (70a) 


Run 89256 
Event 224212 


Multi-electron Event 2 (2e) 

E- P,= 43.9 ± 0.8 GeV p™'*'*=1.9 ± 1.8 GeV 
Mi2 = 130.0 ± 2.6 GeV Alf^ = 129.3 ± 2.4 GeV 


ei 
e2 


132.4 ± 4.3 
82.4 ± 1.8 


28.6 ±0.1 
48.4 ± 0.2 


8.73 ± 0.06 
-171.50 ±0.03 


undetermined 
-(6a) 


Run 254959 
Event 17892 


Multi-electron Event 3 (2e) 

E-P^= 57.3 ± 1.4 GeV P™*««=3.5 ± 2.0 GeV 
Mi2 = 1 12.5 ± 2.4 GeV = 109.5 ± 1.0 GeV 


ei 

62 

fwd em cluster 
photon 


96.9 ± 3.3 
46.1 ± 1.1 

jQ+mOa 

1.1 ±0.1 


34.6 ± 0.3 
80.1 ±0.9 
4.5 ± 0.1 
132.0 ±5.5 


52.66 ± 0.02 
-125.62 ±0.01 
-132.7 ± 1.0 
39.8 ±7.3 


+ (10a) 
+ (15a) 
undetermined 



" this error includes the uncertainty due to energy loss in the beampipe 


Run 168058 
Event 42123 


Multi-electron Event 4 (3e) 

E- P,=55J ± 1.4 GeV p™''^=l.l ± 0.8 GeV 
M12 = 137.4 ± 2.9 GeV M/2* = 138.8 ± 1.2 GeV 


ei 
e2 
es 


35.8 ±0.9 
173.0 ±5.5 
44.8 ± 1.7 


115.6 ±0.9 
6.6 ±0.1 
21.8 ±0.2 


-5.98 ± 0.02 
-159.1 ±0.5 
139.10 ±0.03 


+ (18a) 
undetermined 
- (12a) 


Run 192864 
Event 123614 


Multi-electron Event 5 (3e) 

E- P^= 53.8 ± 1.4 GeV p™*'^=0.7 ± 0.6 GeV 
M12 = 118.1 ± 2.6 GeV = 121.9 ± 0.6 GeV 


ei 
e2 
es 


138.9 ± 4.5 
28.1 ± 0.8 
35.3 ± 1.5 


10.2 ±0.1 
134.7 ± 0.3 
26.6 ± 0.1 


44.1 ±0.6 
-95.85 ± 0.01 
172.71 ±0.05 


undetermined 
+ (25a) 
+ (5a) 


Run 267312 
Event 203075 


Multi-electron Event 6 (3e) 

E- P,= 57.4 ± 1.6 GeV p™*'^=2.4 ± 0.8 GeV 
M12 = 134.7 ± 3.1 GeV M/2" = 132.3 ± 1.4 GeV 


ei 
e2 
es 

hadrons*" 


186.0 ±5.9 
25.5 ± 0.8 
8.5 ± 0.5 
123.2 ± 6.7 


7.11 ±0.05 
148.8 ± 0.2 
69.7 ± 0.3 
6.1 ± 1.1 


-71.3 ± 0.4 
120.25 ± 0.02 
164.90 ± 0.01 

53.5 ± 1.1 


undetermined 
+ (32a) 
+ (57a) 


^ mass of the visible hadronic system: 24.0 ± 2.5 GeV 



Table 4: Reconstructed kinematics of the six multi-electron events with M12 > 100 GeV (see 
text). E is the particle's energy and 6 and its polar and azimuthal angles, respectively . The 
charges of the electron candidates are also given, if they are measured with a significance of 
better than two standard deviations. The six events were recorded in e+p collisions. 
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Electron misidentification study using a DIS event selection 



No track requirement 
• HI data 



a) 



a. 



^3 Pair production 
(GRAPE) 

Compton 




80 100 




(GeV) 



60 

p e-, 



80 100 

Photon conversion study using an elastic Compton event selection 



10 





1 2 3 4 5 
Number of associated tracks 



10 20 30 40 50 60 70 
Track start radius (cm) 



Figure 2: (a and b) Distributions associated with misidentified electrons for a DIS event se- 
lection. Transverse momentum spectrum of second electrons identified with either no track 
requirement (a) or a loose track requirement (b) compared with expectations, (c and d) Study 
of photon conversions using an elastic Compton event selection: number (c) and starting radius 
(d) of the tracks associated with the photon candidates compared with expectations. 
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- • HI Data 115 pb"' 


Pair 




2e 


Production 






(GRAPE) 






-[X3 Dis 






+Compton 



















50 

E-P, (GeV) 




pmiss ^^^y^ 



phadrons ^^.^^^ 




E-P^ (GeV) 



pmiss ^^^^^ 
T 



phadrons ^^^^ 
T 



Figure 3: Distributions of E — P^, p™*'' and p^"-^"^""^^ for events classified as "2e" (top) and 
"3e" (bottom), compared with expectations. 



18 




Figure 4: Electron transverse momentum distributions for events classified as "2e" (top) and 
"3e" (bottom) compared with expectations. Electrons are ordered according to decreasing trans- 
verse momentum, P!p > P^'^\ 
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50 100 150 
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Figure 5: Distribution of the invariant mass M12 of the two highest Pt electrons compared with 
expectations (left) and the correlation of M12 with the scalar sum of the Pt's of the electrons 
(right) for events classified as "2e" (top) and "3e" (bottom). The bold dots on the right plots 
represent the data while the small points represent the pair production (GRAPE) prediction for 
a luminosity 500 times higher than that of the data. 
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100 200 

Mj23 (GeV) 

Figure 6: Invariant mass distributions of electron pairs 1-3 and 2-3 (top left and right) and of 
the tri-electron system (bottom right) for events classified as "3e". 
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HI Data 
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Figure 7: Cross section measurements in a restricted phase space dominated by the photon- 
photon process as a function of the leading electron transverse momentum (P^^), the invariant 
mass of the electron-positron pair (M12) and the hadronic transverse momentum (^p^adronsy 
The differential cross sections are averaged over the intervals shown. The inner error bars on 
the data points represent the statistical error, which dominates the measurement uncertainty. 
The outer error bars show the statistical and systematic uncertainties added in quadrature. The 
bands represent the one standard deviation uncertainty in the SM prediction. 
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Figure 8: Displays of the three "2e" events (left) and the three "3e" events (right) with 
Mi2 > 100 GeV in the R — z view. The reconstructed tracks and the energy deposits in the 
calorimeters are indicated (see text and table|4l). The beam positrons enter the detector from the 
left and the protons from the right. 
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